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Abstract 

The FNAL experiment E773 has measured the phases $+_ = 43.35° ± 0.70° ± 
0.79° and $00 - *+- = 0.67° ± 0.85° ± 1.1° of the CP violating parameters ??+_ 
and ?7oo in the decay of neutral kaons into two charged or neutral pions. These 
preliminary results test CPT symmetry and show no evidence for a violation. In 
addition we present a preliminary measurement of Am = (0.5286±0.0029±0.0022)- 
qI^ 10^° h/sec and ts = (0.8929 ± 0.0014 ± 0.0014) • 10" 

. statistical and the second errors are systematic for the above results. 



^ 



*^ For a long time it was thought that all interactions are invariant under the oper- 

^O ations parity P, time reversal T and charge conjugation C. It was therefore a surprise 

^-^ when in 1957 three experiments [1, 2, 3] found that parity is violated in the weak in- 
teraction. Subsequently in 1964, another experiment found that the combination CP is 

(^ also violated in weak decays of the neutral kaon at the 10"'^ level [4]. 

<^ For the combined operation CPT there is a strong theoretical prejudice stated in the 

L^ CPT theorem [5]: every local field theory (e.g. the standard model) is CPT symmetric. 

■^ To date, no violation of CPT has been reported. 

Q\ Since the strangeness eigenstates K^ and K^ can decay into the same final states 

^ the physical states are linear combinations: 

Dh Ksl = ^ ((1 + e)K° ±(1- e)K°). (1) 

^ V ^ M ; 

CP violation in the neutral kaon decay is parameterized by the measurable ratio of 
violating to conserving amplitude: 

A(i^L^7r+7r-) ^$ 

^+- = A(if5-7r+7r-)='^+-'' ' ^') 

and a similar definition of t/qo for neutral pion decays. If CPT is not violated one can 
show that to within 0.1° [6] 

-1 2Am , , 

$+_ » $e ~ ^sw = tan -^, (3) 



*Talk given at the Conference on the Intersections of Particle and Nuclear Physics, St. Petersburg, 
Florida, June 1994. 
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and Rutgers University 



A$ = $00 - $+- ~ 0. (4) 

with Am = rriL — ms and AT = Ts — T^. 

In E773 a beam of K^ impinged on a regenerator consisting of plastic scintillator. 
After the regenerator the beam was a coherent superposition \Kl,) + p\Ks) with p being 
the amplitude for forward regeneration (scattered kaons were treated as background). 
The decay rate into two pions as a function of proper time after the regenerator is then: 

By fitting the decay distribution for neutral and charged pions the phases $_| and 

A$ can be extracted from the interference term. In addition we can fit for the mass 
difference Am and the Ks life time ts = 1/Ts- 

Experimentally, p follows a power law [7, 8] as a function of kaon energy except for 
a correction due to decays of Ks inside the regenerator (/) ~ p"^ g{p))- If one further 
assumes that the elastic kaon-nucleon scattering amplitudes are analytic functions then 
the phase of p can be related to the power law constant a using dispersion relations. 
Thus the regeneration amplitude can be parameterized by, say, its magnitude at 70 GeV 
and a. 

The experiment E773 took data for three months in the 1991 fixed-target run at 
FNAL. The apparatus (Fig. 1) is essentially the same as for the 1987-88 run (experi- 
ment E731B [8]). Two K^ beams are produced by a 800 GeV jc proton beam striking 
a beryllium target. The detector components are located approximately 120 m down- 
stream of the target. While E731 used one regenerator in one of the beams, E773 ran 
with two regenerators, one at z = 117.14 m (the upstream regenerator, UR) and one 
in the other beam at z = 128.42 m (the downstream regenerator, DR). The UR is 1.2 
and the DR is 0.4 interaction lengths long. Each of them toggled every minute between 
the two beams. Since the regenerators were constructed out of scintillator, i.e. active, 
inelastic interactions of kaons or neutrons in them could be vetoed. 

The decay region extended from the regenerators to the end of the vacuum tank 
at z = 159 m. Inside the vacuum tank and throughout the spectrometer were photon 
vetoes at 11 different z positions to detect charged and neutral particles leaving the 
fiducial region of the detector. Following the vacuum tank was a charged spectrometer 
consisting of four drift chambers and a magnet, which was used to find charged tracks 
and to measure their momenta. The momentum resolution was cjplp = 0.45 (l®p/40)% 
with p in units of GeV jc. After the spectrometer there were two scintillator planes 
with horizontal or vertical segmentation. The trigger for the two-body charged decays 
K -^ 7r"'"7r~ required up-down and left-right hits in the hodoscopes and drift chambers. 
For part of the run there was another trigger hodoscope at z = 141 m which was 
removed, after a trigger upgrade, to extend the decay region for charged decays. 

Following the hodoscope was a circular array of 804 lead glass blocks (each block 
5.8 cm by 5.8 cm and 18.7 Xq deep) with holes in the middle of the array for the 
two beams. This was used as an electromagnetic calorimeter for the neutral decays 
K -^ 7r°7r°. The trigger required at least 25 GeV in the calorimeter and 4 or 6 clusters 
found by a hardware cluster finder. 

All electromagnetic shower energy was contained in the lead glass and the lead wall 
(21.4 Xq) downstream of it. A scintillator plane further downstream was therefore used 
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Figure 1: Schematic detector picture of experiment FNAL E773. 



to veto hadronic showers. Downstream of this plane were 3 m of iron and another 
scintillator plane. Muons from K -^ tt/xz/ decays would penetrate this material and 
could thus be detected and vetoed in the charged mode trigger. 

Drift chamber calibration constants (transverse position and rotation of chambers 
and time-to-distance relations) were generated about once a day by using tracks from 
special runs with the magnet off and "two-track" decays. For the lead glass calibration 
a large sample of electrons from K -^ 'Kev decays was used. Because of increasing 
radiation damage throughout the run the opacity of the glass changed. Blocks that were 
hit most were therefore calibrated about once a day. The gain of the photo tube/ADC 
and the light absorption coefficient for every block was determined. 

The lead glass was calibrated using electrons. Because of absorption the light output 
of a photon shower at a given energy was different from that of an electron shower. To 
correct for the size of the electron-photon cluster differences we relied on a simulation 
which was fine tuned by comparing the Monte Carlo and data z distributions oi K ^> 
7r°7r° and K -^ 37r° around the regenerators. The z position of this sharp edge was 
sensitive to small changes in the energy scale and a small energy dependent correction 
(about 0.3%) was applied. 

To reconstruct K -^ tt+tt" events, two tracks were found in the x and y view and 
then matched to clusters of energy in the lead glass. The ratio of cluster energy to spec- 
trometer momentum E/p had to be less than 0.8 to reject electrons from semileptonic 
decays K -^ 'Kev. Track and vertex quality cuts as well as aperture cuts were applied. 
A minimum momentum requirement ensures that a muon would have penetrated the 
muon filter and fired the muon veto. Events with a pi: mass within 6 MeV jc^ of the A 
mass were rejected to reduce A ^ pTr background. K -^ 7r"'"7r~7r° events were rejected 
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Figure 2: Ring number distributions oi K ^ tt tt for the downstream (part a) and 
upstream regenerator (part b). Shown are the data (histogram) and a Monte Carlo 
simulation of the background from scattered kaons( circles). The cut was at 120. 



by requiring the 7r"'"7r~ mass to be close to the kaon mass. Finally, kaons that scattered 
in the regenerator with a transverse momentum pj^ > 250 (^MeV/cy were discarded. 
After all cuts the background was quite small (0.3% and 0.8% for the UR and DR beam, 
respectively). 

For K -^ 7r°7r° decays the four clusters in the lead glass were combined into two 
pairs. Assuming two 7r° decays the z decay vertices of the pions were calculated. With 
four photons there are three combinations. If the two z values agreed with a x^ < 4 
the event was kept. Further requirements were: four cluster mass between 474 MeV/c'^ 
and 522 MeVjc^, off-line cut on photon veto energies, and cut on the shape of the 
electromagnetic cluster (to ensure that each cluster came from a single photon). All 
three cuts reduced the K -^ 7r°7r°7r° background. Events with out-of-time energy due 
to accidental activity were rejected. The four photon center of energy was the position 
where the kaon would have hit the calorimeter if it had not decayed. The maximum 
distance between this position and the center of the beam in the x or y view (called 
"ring number") had to be small. For a kaon that scattered in a regenerator it was 
likely that its position at the glass was outside the beam profile and would therefore be 
rejected. 

The background to the neutral mode came from three different sources. For K^ -^ 
7r°7r°7r° with two photons escaping detection or fusing, a Monte Carlo simulation was 
used. The normalization was taken from the mass side bands since all events in the 
data that reconstruct outside the kaon mass and do not come from beam interactions, 
i.e. do not reconstruct at the z positions of vacuum windows, were K -^ 37r° decays. 
The background fractions were 0.6%(1.4%) for the UR(DR) beam. 



c 400000 — 

a3 350000 - 
Cl : 

-^ 300000 7 
CD ; 

L^ 250000 7 

200000 

150000 7 

100000 7 

50000 7 



O 



o 
o 



7^ 

1.06 7 
1,04 7 

1.02 7 

:+ 



a) 




+ i, A 



1 r+%#+^^^+++V+++V 



0.98 7 
0.96 7 
0.94 - 



120 130 140 150 
Z from target in m 



28000 



24000 



Q_ 

-I-' 

C 

^ 20000 



O 



o 

q 
o 



16000 

12000 

8000 

4000 


1.06 

1.04 

1.02 

1 

0.98 

0.96 

0.94 



b) 




120 130 140 150 
Z from target in m 



Figure 3: Decay distribution of the upstream regenerator beam (a) K^ -^ T^ev and (b) 
Kl -^ 37r° decays. The upper plots show data (histogram) and Monte Carlo (circles) 
while the lower plots show the ratio of data over Monte Carlo. 



The background from beam interactions in the vacuum windows or air around the 
regenerators and the regenerators itself was subtracted by linear interpolation in the 
four photon mass distribution. The size of this background was about 0.1%. 

The third class of background was due to scattered kaons. To estimate the amount 
of background that was left after the ring number cut a Monte Carlo simulation of 
the scattering was used. The input to the simulation was taken from the K -^ 7r"'"7r~ 
decays, since the kaon pj^ could be calculated for these events. The kaon could either 
scatter inelastically which results in a rather flat p^ distribution (exp( — 5p|i) with p^ 
in units of GeV'^/c^) or diffractively. The diffractive K^ -^ Kl scattering was param- 
eterized by exp( — 58p|i) while the diffractive Ki, -^ Ks scattering was parameterized 
by exp( — 90p|i). The exponents were the same for both regenerators. Together with 
two normalizations, one for each regenerator, the diffractive background was thus pa- 
rameterized by four constants in the Monte Carlo which predicted the z, energy and 
Pj^ distributions for all decay modes quite successfully. Fig. 2 shows the ring number 



distribution for 7r°7r° decays and the Monte Carlo simulation of the background. The 
background fractions were 3.0% for the upstream (thicker) and 9.7% for the downstream 
(thinner) regenerator. 

Understanding the acceptance of the detector was crucial for this experiment since 
the reported measurements are obtained by fitting the acceptance corrected decay dis- 
tributions to the expected rate (5). The acceptance was checked by comparing high 
statistics data and Monte Carlo z distributions of Ki, decays into 37r° and Trez/ for the 
7r°7r° and tt+tt" decay mode, respectively. These decay modes were not sensitive to the 
simulation of the regeneration. Fig. 3 shows that the acceptance is understood to within 
1%. 

Fig. 4 shows as an example the corrected distribution of if ^ 7r"'"7r~ for one kaon 
energy bin. Three different kinds of fits to the background subtracted and acceptance 
corrected decay distributions were performed. In all fits the two quantities parameteriz- 
ing the regeneration amplitude p (|/o(70 GeV')\ and a), normalizations for the kaon flux 

and two shape parameters for its energy dependence were floated. In the first fit $_| 

was floated in addition. In the fit for Am and ts-, $h was fixed to ^^^ (assuming CPT 

invariance, eqn. 3). For A$, both phases, |e'/e| (3e' = 7y_| — t/qo)) a normalization for 

the relative number of 7r°7r° to 7r"'"7r~ events and two more energy dependend parameters 
for the K -^ 7r°7r° flux -^ were floated. 

The fit result for $_| and its dependence on Am and ts (for one Particle Data 

Group (PDG) standard deviation [9]) is: 



$+_ = 43.35° ± 0.70°(5iai) ± ^.Td\syst) 

„Ts- 0.8922 •10-i°sec „ Am - 0.5286 • 10i°;i/sec ,, 

0.0020 -lO-iOsec ' 0.0024-1010^560 '^' 



Note that the central value for ts is taken from PDG but that the Am value is our best 
fit which is motivated by the argument described below. 

The systematic errors for the fits are listed in table 1. For the A$ fit the systematic 
error due to the minimal cluster energy cut as mentioned below is part of the systematic 
due to the nonlinear lead glass response. 

In the Ts and Am fit we get: 

Am = (0.5286 ± 0.0029(siai) ± 0.0022(s2/si)) • 10^°;i/sec, (7) 

TS = (0.8929 ± 0.0014(siai) ± 0.0014{syst)) ■ 10"^°sec. (8) 

The result for Am is very interesting since it deviates from the PDG value AmpDc = 
(0.5351 ± 0.0024) • 10^^%/ sec [9], but confirms the most recent measurement from E731 

[13]. This is relevant for the CPT test $_| ~ ^sw since the value of ^sw depends on 

Am and because there is a strong correlation between Am and $_| in the phase fit. 

This is true for all published measurements. The current PDG values are $_| ^pdg = 

46.6° ± 1.2° and ^swPDG = 43.73° ± 0.15° which hints at CPT violation. If however 



Because of acceptance effects not simulated in the Monte Carlo there was a small mismatch in the 
energy spectrum for tt tt and tt 7r~ decays. 



LO 

d 

CD 
CL 



C 

> 



10' 



10' 



10- 



Acceptance Corrected Data 
Prediction with Interference 
Prediction w/o Interference 




1..^ 



yj^t^ 



120 125 130 135 140 145 150 155 

Z decay In m from target 



Figure 4: Background subtracted and acceptance corrected decay distribution of if ^ 
7r"'"7r~ for the upstream regenerator beam in the energy bin 40 — 50 GeV (crosses) and 
the predictions from the fit with and without the interference term. 
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Figure 5: Previous measurements of $_| and this result (open square). The open circles 

are the results used to determine $_| ,pdg, the solid points are the same measurements 

corrected using our Am and the solid squares are more recent results. The band is 
the ±1 (7 band of ^sw The corrected values all agree with each other and with ^sw 

Thus the test $_| ~ ^sw shows no hint of CPT violation anymore. The ts and Am 

dependences are included in the errors. The references to the results are in time order 
[10, 11, 12, 13, 14]. 
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Table 1: Systematic error estimates for the fits. The individual sources were added in 
quadrature to give the total error. 





$+_ fit 


Am/rs fit 


A$fit 


in degree 


in W^hlsec 


in 10-^"sec 


in degree 


regenerator positions, lengths 


0.1 


0.0008 


0.0004 


- 


charged mode acceptance 


0.57 


0.0017 


0.0010 


0.2 


neutral mode acceptance 


- 


- 


- 


0.2 


background subtraction 


0.2 


0.0005 


0.0002 


0.6 


deviations from the power law 


0.5 


0.0011 


0.0008 


- 


for the reg. amplitude 










lead glass energy scale 


- 


- 


- 


0.6 


lead glass energy resolution 


- 


- 


- 


0.5 


lead glass nonlinear response 


- 


- 


- 


0.5 


total 


0.79 


0.0022 


0.0014 


1.1 



previous experimental results are corrected for their Am dependence to our lower value 
of Am the discrepancy vanishes. This is shown in Fig. 5. The open circles are the results 

used to determine $_| ^pdg while the solid points are the corrected measurements and 

more recent results (solid squares). The hatched band shows ^sw using our central 
value for Am. Note that after the correction there is perfect agreement between ^sw 
and $_| and among all experiments. 

The preliminary Am result from CPLEAR, presented at this conference [15], is 
Am = (0.5318 ± 0.0044(siai) ± 0.0010(s2/si) ± 0.0030(MC)) • 10^° hisec, consistent with 
our and E731 measurements. 

For the A$ fit we used 70% of the data for the preliminary analysis. The fit result 
was sensitive to the cut on the minimal cluster energy in the lead glass. Changing the 
cut from 2.2 GeV to 4.0 GeV changed the result by 0.9° which is not yet understood. 
For the preliminary result we have decided to take the average of the two fits as the 
central value and to assign 0.5° as a systematic error due to this effect. The preliminary 
result is: 



A$ = 0.67° ± 0.85°(siai) ± l.l°[syst). 



(9) 



We find that A$ is consistent with zero (eqn. 4) and our values of ^sw and $_| agree 

(eqn. 3). Therefore both tests find that CPT symmetry is valid to the accuracy of our 
experiment. 

In addition to the CPT tests, E773 has the best measurements of the CP violation 
parameter in the radiative decay K -^ 7r"'"7r~7. Our preliminary results are: 



l^+- 



(2.414 ± 0.065(siai) ± 0M2{syst)) ■ 10"^, 



$ 



+ -7 



45.5° ± 3.6°{stat) ± 2A°{syst). 



(10) 
(11) 



Thus 7y_| y for the radiative decay (inner Bremsstrahlung) is consistent with being the 

same as 7y_| [9]. These measurements establish an improvement in precision by more 

than a factor of 3 relative to the previous results [16]. 



The results presented comprise the thesis work of Roy A. Briere (University of 
Chicago), John N. Matthews (Rutgers University) and Bernhard Schwingenheuer (Uni- 
versity of Chicago). 
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